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Abstract 

A  semi-empirical  model  of  helical 
magnetocumulative  generators  (MCGs)  has  been 
developed  and  will  be  presented  in  this  paper.  A 
single  adjustable  variable,  called  the 
characteristic  time,  has  been  identified  and  its 
physical  significance  of  this  variable  will  be 
discussed.  The  model  has  been  applied  to  four 
different  helical  MCGs  and  the  calculated  results 
are  in  good  agreement  with  the  experimental 
results  (See  Paper  P1-E29  of  this  conference). 

I.  INTRODUCTION 

The  model  presented  in  this  paper 
describes  the  time  dependent  changes  in  the 
output  current  of  helical  generators.  This  semi- 
empirical  model  appears  to  be  general  enough  to 
serve  as  an  initial  design  code.  Its  simplicity  and 
speed  make  it  practical  for  doing  parametric 
studies  of  various  types  of  helical  generators 
with  different  types  of  loads. 

II.  HELICAL  MCG  MODEL 

All  the  above-mentioned  generators  can 
be  represented  by  a  series  LR-circuit  described 
by  the  differential  equation: 


^■[(4+l1);]+«/= o,  (i) 

where  Lc  and  Ll  are  the  inductance  of  the 
generator  and  load,  respectively,  R  the  resistance, 
I  the  current,  and  t  is  time.  The  inductance  of  the 
generator  and  the  resistance  will  vary  in  time, 
while  the  inductance  of  the  load  remains 
constant.  Integrating  this  equation  yields: 

I  =  I0e-g(,\  (3) 

where 


*w=J 
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and  I0  is  the  initial  current. 

The  inductance  of  the  coil  in  the  MCG 
is  determined  by  using: 


/yrA2 

/ 


(5) 


where  [i0  is  the  permeability  of  free  space,  N  is 
the  number  of  turns,  l  is  the  length  of  the  coil, 
and  rc  and  ra  are  the  radii  of  the  coil  and  the 
armature,  respectively,  and  the  number  of  turns 
in  the  coil  can  be  written  in  terms  of  its  length 
and  pitch,  N  =  Up.  As  the  armature  expands  and 
successively  short-circuits  the  turns  in  the  coil, 
the  length  of  the  coil  is  shortened.  Assuming 
that  this  can  be  described  by  the  equation: 
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where  x  is  a  characteristic  time,  and  integrating 
Eq.  7  yields: 


i 


l  =  l0e  T . 
The  inductance  is  defined  by: 


P 

The  resistance  can  be  written  in  the 


(8) 

(9) 


form: 


R  =  — ,  (10) 

oA 

where  a  is  the  electrical  conductivity,  AR  is  the 
area  of  the  region  contributing  to  the  resistance, 
and  /  is  the  length  of  the  coil.  Substituting  Eq.  8 
into  this  equation  yields: 

l  ~- 

R  =  —9—eT=RQeT.  (11) 

°Ar 

Substituting  Eqs.  9  and  1 1  into  Eqs.  3 
and  4  yields: 


/  =  /o  exp 
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It  is  assumed  that. 


RqT  . 

— —  <&l,  (13) 

A> 

for  the  devices  under  consideration.  This 
condition  indicates  that  initially  the  current  in  the 
MCG  is  increasing  as  predicted  by  Eq.  12  and 
verified  experimentally. 

From  a  physical  point  of  view,  it  is  clear 
that  at  some  time  to,  the  current  must  stop 
increasing  and  start  to  decrease.  The  time  at 
which  this  occurs  should  be  related  to  the  time  at 
which  the  detonation  ends.  This  means  that  Eq. 
12  must  be  modified.  Empirically,  it  can  be 
shown  that  the  experimental  data  is  represented, 
to  a  reasonable  degree  of  accuracy,  by  the 
equation: 
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where  t  >  t0  and  b  =  1.8  x  105  sec'1.  It  should  be 
noted  that  when  t  =  t0,  Eq.  14  reduces  to  Eq.  12. 

For  the  purposes  of  this  study,  it  has 
been  determined  empirically  that  the  initial 
resistance  for  each  of  the  generators  is  given  by 
the  relation: 


Ro  =  ~ tA  U5) 

A 

where  la  is  the  length  of  the  armature.  The 
constant  in  Eq.  15  has  the  units  of  electric  field 
(V/m)  and  911  a  is  an  equivalent  voltage. 

The  time  t0  is  given  by  the  relation: 

la 

t0=—,  (16) 

vD 

where  vD  is  the  detonation  velocity.  In  other 
words,  to  is  the  time  at  which  the  detonation 
terminates. 

The  initial  inductance  is  given  by  a 
series  of  equations  of  the  same  form  as  in  Eq.  9, 
that  is: 


L0=JuMrc2~ra2)Y\,  (17) 


*=i  Pk 

where  K  is  the  number  of  induction  coils. 


III.  Characteristic  Time 

In  the  above  analysis,  the  only 
parameter  adjusted  empirically  in  order  to  fit  the 
calculated  data  to  the  experimental  data  is  the 
characteristic  time,  x.  The  values  of  the 
characteristic  time  for  the  helical  generators 
considered  are:  Mark  IX  -  48.0  ps,  EF-3  -  31.0 
ps,  Flexy  I  -  29.2  ps,  and  Ranchito  -  3 1  ps.  The 
calculated  results  are  very  sensitive  to  the  value 
of  x.  One  of  the  objectives  of  this  paper  is  to 
begin  to  understand  what  physical  factors 
influence  this  parameter.  Some  of  these  factors 
include  geometrical  sizes  and  shapes,  material 
properties,  and  detonation  characteristics. 

As  a  first  step  in  determining  what 
factors  influence  x,  it  is  written  in  the  following 
form: 


(18) 


where  lc  and  vc  are  a  characteristic  length  and 
velocity,  respectively.  As  the  detonation  wave 
proceeds  down  the  armature,  the  armature 
expands  with  some  velocity  until  it  hits  the  coil. 
Thus,  the  vertical  distance  that  the  armature 
travels  is  lc  =  rc  —  ra.  The  expanded  armature 
forms  an  angle  a  with  the  axis  of  the  generator. 
Therefore,  the  characteristic  velocity  is: 

vc  =  fvD  tana,  (19) 

where  the  fraction /has  been  introduced  to  take 
into  account  that  only  a  fraction  of  the  detonation 
energy  is  transferred  into  motion  of  the  armature. 
Therefore,  the  characteristic  time  is: 

r  —r 

T  =  — ^ .  (20) 

fvD  tan  a 

The  value  of/ depends  on  the  mass  of 
the  liner  and  the  mass  of  the  explosives  available 
to  provide  the  energy  for  bending  and  moving 
the  armature.  As  a  first  approximation,  the 
following  relation  will  be  used: 

'  c  y 

/=  - M—Pr  ,  (2D 

1  +  0.5  — 
l  M 

where  C  and  M  are  the  mass  of  the  explosive  and 
the  armature,  respectively.  This  relationship  was 
developed  by  Gurney  [1]  in  his  analysis  of  the 
velocity  of  fragments  from  an  explosion. 

In  the  Mark  IX  experiments,  60  kg  of 
PBX  9501  was  used  and  the  armature  was  made 
of  copper.  Using  the  dimensions  of  the  armature 
given  in  [2],  the  calculated  mass  of  the  armature 
is  71.9  kg,  which  means  that  C/M  =  0.834  and/= 
0.767.  Using  these  values  and  a  characteristic 
time  of  48  liscc,  the  angle  a  is  calculated  to  be 
15.7°.  The  measured  value  given  in  [2]  is  14°. 

For  the  Flexy  I,  the  armature  was  made 
from  aluminum.  Using  the  dimensions  given  in 
[3],  its  mass  is  calculated  to  be  8.29  kg.  The 
mass  of  the  explosives  was  given  as  15  kg. 
Therefore,  C/M  =1.81  and/=  0.975.  For  a 
characteristic  time  of  29.2  liscc,  the  angle  a  is 
12.8°.  The  measured  value  given  in  [3]  is  12°. 

For  the  EF-3,  the  mass  of  the  explosives 
used  was  not  given,  but  the  angle  a  was  given  in 
[3]  as  13°.  The  armature  was  made  of  copper 
and  its  mass  is  calculated  to  be  9.75  kg. 
Assuming  that  the  mass  of  the  explosives  is  7  kg, 
it  is  found  that  C/M  =  0.72  and/=  0.727.  Using 
these  values  and  a  =  13°,  the  characteristic  time 
is  calculated  to  be  3 1  liscc. 


IV.  SUMMARY 

A  semi-empirical  model  of 
helical  MCGs  has  been  developed  and 
benchmarked  by  using  it  to  calculate  the 
output  parameters  of  generators  having  a 
wide  range  of  design  and  operating 
conditions.  In  deriving  the  model,  a 
sensitive  adjustable  parameter,  called  the 
characteristic  time,  was  identified.  This 
parameter  is  influenced  by  factors  such 
as  geometrical  sizes  and  shapes,  material 
properties,  and  detonation 
characteristics.  Relations  were 
developed  that  show  how  these  factors 
affect  the  characteristic  time  and,  thus, 
operation  of  the  generators. 
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